Surface Phase Transitions Induced by Electron Mediated Adatom-Adatom Interaction 
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We propose that the indirect adatom-adatom interaction mediated by the conduction electrons of a metallic 
surface is responsible for the V3 x V3 « 3 x 3 structural phase transitions observed in Sn/Ge (111) and 
Pb/Ge (111). When the indirect interaction overwhelms the local stress field imposed by the substrate registry, 
the system suffers a phonon instability, resulting in a structural phase transition in the adlayer. Our theory is 
capable of explaining all the salient features of the V3 x V3 <=> 3 x 3 transitions observed in Sn/Ge (111) and 
Pb/Ge (111), and is in principle applicable to a wide class of systems whose surfaces are metallic before the 
transition. 
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PACS numbers: 68.35.Bs,68.35.Rh,73.20.At,71.45.Lr 



Over the years, a great deal of efforts have been devoted 
to experimental studies of structural phase transitions at sur- 
faces. One compelling example is the V3x V3 <=> 3x3 transi- 
tions in the 1/3 monolayer of Pb or Sn on Ge( 1 1 1 ) directly ob- 
served by the Scanning Tunneling Microscopy (STM) 
(see Fig. These studies have stimulated an active line 
of theoretical research, yet to date, the precise nature and 
the underlying mechanism of such transitions are still highly 
controversial. The original paper attributed this transition to 
a Charge Density Wave (CDW) driven by two-dimensional 
Fermi surface nesting Ql @|- Subsequent papers have at- 
tributed the transition to a Kohn Anomaly |4fl, bond den- 
sity waves a pseudo Jahn-Teller transition |6, 7], a sur- 
face Mott insulator (H, dynamical fluctuations ff, a soft 
phonon 1 10] and most recently to disproportionation II ill . Al- 
though all these theories are successful to certain extent, a uni- 
fied microscopic picture is yet to emerge. The issue becomes 
even more intriguing after the observation of the delicate role 
of defects in the transition iflil [rUl . 

In this Letter, we present a new mechanism for sur- 
face phase transitions, which places central emphasis on 
the indirect adatom-adatom interaction mediated by the two- 
dimensional conduction electrons of a metallic surface. In 




J3XJ3 3x3 
) Up Sn 9 Down Sn 3 1st Layer Ge 9 2nd Layer Ge 

Figure 1 : Structures of V3 x V3 (left) and 3x3 (right) phases of 
Sn/Ge(lll). In the experiment, one third monoatomic layer of tin 
was deposited onto a Ge(lll) surface, and a gradual transition from 
the room-temperature flat V3x V3 phase to the low-temperature 3x3 
phase was observed [lj|2jj. In the new phase, one Sn adatom moves 
up and two move down in each 3x3 unit cell. 



this theory, when the conduction-electron mediated adatom- 
adatom interaction overwhelms the local stress field imposed 
by the substrate registry, the system suffers a phonon insta- 
bility, leading to a structural phase transition. The theory is 
capable of explaining all the salient features of the transition 
observed in Sn/Ge (111), including the appearance of surface 
charge density waves and the delicate role of surface defects. 
It also predicts the existence of a glassy phase in such systems. 

The electron-mediated adatom-adatom interaction origi- 
nates from the tendency of the conduction electrons to screen 
external disturbances. In the case of adlayer systems such as 
Sn/Ge (111), the dangling bond electrons of the adatoms form 
a quasi two-dimensional (2D) electron gas at the surface. The 
movement of an adatom disturbs the electron gas, which re- 
sponses in the form of Friedel oscillations in its density. Such 
charge corrugations propagate at the surface to reach other 
adatoms, thereby establishing an indirect interaction between 
adatoms. We note that the basic physical concept empha- 
sized here has found important applications in several other 
fields, such as the RKKY interaction between spins 1 14| and 
chemisorption of adatoms at metal surfaces 1 15Ul6l 

We stress that, the central ingredients of the picture de- 
veloped here are the existence of surface conduction elec- 
trons and their capability of coupling with the displacement of 
the adatoms (electron-phonon coupling). Both are evidently 
present in Sn/Ge (11 1) and Pb/Ge (111). As shown in Fig.[0 
each Sn adatom bonds to three Ge substrate atoms directly 
under it, leaving one bond dangling. The electrons in the dan- 
gling bonds are localized in the surface and form a 2D elec- 
tron system with a narrow (~ 0.5 eV) half-filled band. When 
an adatom is displaced, the angles between the dangling bond 
and its three saturated Sn-Ge bonds have to adjust accordingly, 
inducing a variation in the s-p hybridization of the electron 
states of the Sn adatom, and a corresponding change in the 
energy of the dangling state Hill . This process can be de- 
scribed by the Holstein electron-phonon coupling 1 19], where 
the dangling state energy of an adatom depends on its dis- 
placement as ed = eo — 0Z, with z being the displacement of the 
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adatom along the direction perpendicular to the surface, and /3 
a positive constant, reflecting that a downward displacement 
induces a higher dangling state energy. 

The total Hamiltonian of the system is written as, 



H = J](ft-Af)4 



t J] C c 



js+ 2 



(1) 



where z, is the displacement of the adatom located at site i, 
and c! (c«) is the creation (annihilation) operator of a dan- 
gling bond electron at site i and with spin s. The hopping 
constant t is between nearest neighbors and can be regarded 
as independent of the adatom displacement because the lead- 
ing correction term is second order in zu The term involving a 
is the elastic energy of the lattice distortion caused by the dis- 
placements of the adatoms, which represents the local stress 
field imposed by the substrate. In Eq.[2 the higher order terms 
in Zi are ignored for the moment. As it will become clear in 
later discussions, Eq.[Qis sufficient for determining the struc- 
tural stability of the system, while the higher order terms are 
important only in stabilizing the system after the system be- 
comes unstable. 

The stability of the V3 x V3 phase can be examined within 
the perturbation theory by expanding the total energy change 
up to the second order in the adatom displacements, 

AE = \ Z ~P {n) Z Zi - \ Z J yfc " z > )2 + 0(z3) • (2) 

> i u 

where («) is the average number of the dangling electrons for 
each adatom, and 7y are the coupling coefficients for the in- 
direct adatom-adatom interaction. In the continuum limit, we 
have, 



8 ^ 2 ZZ 

8n 2 yS : 



exp[-i(q-k).(Ri-Rj)] 



F(k F \Ri - Rj\) 



(3) 



where e q {e^) is the energy dispersion of the electron band, 
€p is the Fermi energy relative to the band bottom and con- 
tains implicitly the hopping constant t. The coefficients 
Jij are the same as for the 2D RKKY coupling 12011 . and 
are oscillatory spatially with an asymptotic dependence of 
- sin(2kfr)/r 2 , as shown in Fig. |2ja). When applied to the 
case of V3 x V3 structure of Sn/Ge(lll), the coupling be- 
tween the nearest neighbors is positive; therefore, the interac- 
tion is anti-ferromagnetic-like, namely, two nearest neighbor- 
ing adatoms tend to displace in opposite directions. 

Keeping only the interacting terms between nearest neigh- 
bors, we can reduce Eq.|5]to 



AE = Ua - 67,) £ ~zf + y J] ziZj + 0(z 3 ) , 



(4) 



(U) 



where zi — Zi — f}/(a(n)). This equation shares the same form 
as the phenomenological charge compensate model (CCM) 




Figure 2: (a) F(x) in Eq.|5|(as a comparison, the triangles are the 
results of tight-binding calculations up to the fifth nearest neighbors 
for the V3x V3 structure). In this system, the Fermi momentum kf » 
^J2np = V4~7r/(3 1/4 a), where p is the density of the surface electrons, 
a is the lattice constant, (b) Phonon dispersion [Eq. |4) along the 
direction TM shown in the inset, with ajj\ = 7. The shaded area 
indicates the unstable phonon modes, of which the most unstable K' 
mode defines the 3x3 periodicity. 



proposed by Melechko et al. l2lll . which has been shown to 
be capable of interpreting most of the experimental STM im- 
ages. Our theory thus provides the microscopic mechanism 
behind the success of the CCM model. We note that the long- 
range terms in Eq.[2]should be included in future investigation 
of large-scale phenomena such as domain wall formation, but 
they do not alter the nature of the phase transitions, which is 
the focus of the present study. 

From Eq. |4] it is evident that the system becomes unstable 
against adatom displacements when the elastic restoring force 
is weak. This can be illustrated with the phonon dispersion of 
Eq.Slll: 



OJ' 



— a - 6 J i + J i ^ cos(k • R/), 



(5) 



where the sum is over six nearest neighbors. When ajJ\ < 9, 
oj becomes imaginary for certain k vectors (Fig. |2>), indi- 
cating that the corresponding phonon modes become unsta- 
ble (phonon catastrophe) and the system is driven into new 
phases. Note that the soft phonon picture suggested in RefUfll 
is a very special case of our model, occurring only when 
a/Ji is exactly equal to 9. Furthermore, the phonon insta- 
bility established here does not require any special properties 
of the Fermi surface such as the Fermi surface nesting, nor 
does it rely on the inclusion of electron-electron correlations, 
which is the focus of most of the earlier efforts with little suc- 
cess mm. 
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(a) a/Ji > 9 (b) 6 < a/Ji < 9 (c) a/Ji < 6 



Figure 3: Schematic images of the stable configurations for different 
parameter regimes, (a) V3 x V3 phase; (b) 3 x 3 phase; (c) kinked- 
line phase. 

Higher order terms become important when the system en- 
ters the unstable regime. On a microscopic level, those terms 
may arise from a variety of sources, such as anharmonicity 
in the elastic energy, finite width of the electron band, higher 
order terms in expansions of the site energy and the hopping 
constant t. The precise form of these higher order terms may 
affect the quantitative details of the phase transition, such as 
the amplitudes of the adatom displacements and preference 
among the nearly degenerate ground states, but the occurrence 
and nature of the transition is totally determined by the intrin- 
sic instability represented in Eq. 0] For this reason and also 
for simplicity, we assume that the higher order terms mainly 
come from the anharmonic terms of the on-site elastic energy, 

i i 

where 5 > and y > 0. S is positive, because displacing an 
adatom out of the surface weakens the bond strengths between 
the adatom and the substrate atoms. 

By adding those higher order terms into Eq. [4] we can de- 
termine the stable configurations. The numerical results are 
shown in Fig. [3] The ratio a/J\ is a crucial deterministic pa- 
rameter for the phase transition. When a/Ji > 9, the ground 
state is the flat V3 x V3 phase. When 6 < a/Ji < 9, the 
ground state is the 3x3 phase in the one up-two down con- 
figuration: in each unit cell of the 3x3 lattice, one Sn adatom 
moves up (zf), while the other two move down (zj,), and the 
displacements satisfy Z| * 2|zjJ . The nearly degenerate one 
down-two up configuration has a higher energy due to the cu- 
bic term in Eq.|6] Such 3x3 patterns were also observed ex- 
perimentally 171 12311 . Note that the 3x3 structure corresponds 
to the most unstable phonon mode in Eq. [5] as demonstrated 
in Fig. 0b). 

When a/Ji < 6, the system shows kinked-line configura- 
tion and behaves like a glass 1 24] . Starting with random initial 
configurations, we always end up with a metastable disordered 
structure as typified in Fig. |3jc), instead of the true ground 
state 3x3. We expect a similar behavior in a realistic exper- 
iment with finite cooling rates: the system is trapped in one 
of the meta-stable states, showing a disordered configuration. 
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Figure 4: Phase diagram of the system. The boundary between the 
glassy phase and the 3x3 ground state is not well defined as indicated 
with the dashed line. T = J\ly, (5 = 0. 

Such behavior is directly related to the property of the single 
body potential of an adatom , 

l -( a -6J l )z^- 6 -z] + 7 -t, (7) 

which develops a double-well shape when ajj\ < 6. As a 
result, the present model can be mapped onto an Ising an- 
tiferromagnet by considering the displacement z; taking dis- 
crete values at the local minima of the double-well potential. 
Such an Ising antiferromagnet on a triangular lattice is known 
to have an exponentially large number of degenerate ground 
states l2*3l . Here, the vibrations of z; around the local minima 
break the degeneracy, leading to the exponentially large num- 
ber of meta-stable states in addition to its true ground state 
3x3. 

The dependence of the stable configurations on ajj\ quali- 
tatively explains why the V3 x V3 o 3 x 3 transition is ob- 
served in the Sn/Ge and Pb/Ge systems, but not in Sn/Si: the 
stronger bond strength between a Sn adatom and a substrate 
Si atom places the system into the regime where a/J\ > 9. 
This is to be verified by future first-principles calculations. 

The structural phase transition manifests itself by an ac- 
companying charge density wave, as observed experimen- 
tally |01 01 • Following the Hellmann-Feynman theorem, we 
obtain the total force acting on an adatom, 

F, = az ; -/?<«,> +0(z 3 ). (8) 

Setting F, = in equilibrium, we have («,■} = («) + (a/fi)z.i + 
0(z 3 ), showing that adatoms displaced upward gain electrons 
while those downward lose electrons. 

We now study the finite temperature behavior of the system 
by applying the Monte-Carlo algorithm to the effective clas- 
sical Hamiltonian shown in Eq.|4] along with Eq.|6] We find 
that the transition between the v3 x V3 phase and the 3x3 
phase is of second-order with a sharply defined boundary, as 
shown in Fig.|4] There exists another boundary between the 
glassy phase and the 3x3 phase, which is the result of the 
glassy states at T — when a/J\ < 6. As is typical in a 
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systems whose surfaces are metallic before the transition. 
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Figure 5: Temperature dependence of the order parameter for the 
defect-free system and the system with 1% defect. The thick solid 
line shows the function A{l-T/T c ) [/2 with T c /T = 0.18, T = Jf/y. 
The parameters are ajJ\ = 8.5, J\ = 1, y = 1.0, and S = 0.1, with a 
simulation system size of 99 x 99. 

glassy system, the boundary of the glassy phase is not well 
defined. The detailed behavior of the V3 x V3 o 3x3 phase 
transition along the temperature axis is shown in Fig. [5] The 
order parameter is chosen as the mean square corrugation of 
the thermo-average positions of adatoms |21]. The defect-free 
system clearly shows a second-order transition from the 3x3 
phase to the V3 x V3 phase as seen in Fig. [5] where the tem- 
perature dependence of the order parameter is well fitted by 
A(l - T/T c ) l t 2 . This behavior is different from the prediction 
of the Ginzburg-Landau type theories such as the CCM, which 
concludes that the critical behavior follows |1 - T/T c \ lElll . 

Figure|5]also shows that the sharp phase transition in a pure 
system is blurred by the presence of a very low concentra- 
tion of substitutional defects. In our calculations, the defect is 
simulated by introducing a constant displacement, 

Az = -Ae/j8. (9) 

where Ae is the site energy difference between a Ge substi- 
tutional defect and a Sn adatom |4]. As seen in Fig.|5| with 
only 1 % defects, the sharp transition in the order parameter 
converts into a crossover behavior, explaining the experimen- 
tal observation that the transition is gradual. Just like in the 
experiments our simulations show that each defect 

induces a local 3x3 patch above T c , and the size of each in- 
dividual 3x3 patch decreases gradually with the temperature. 

In conclusion, we have presented a complete theory to un- 
derstand the surface phase transitions observed in Sn/Ge (111) 
and Pb/Ge (111). Although we focus on these specific sys- 
tems, the central ingredients of the theory, namely, the elec- 
tron mediated indirect interaction and the resulting phonon 
instability, are conceptually much more general. These ba- 
sic ideas, once properly applied, will enable us to gain better 
understanding of structural phase transitions in a wide class of 
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